Pseudogap above Tc in a model with dx2_y2 pairing 
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We study the anomalous normal state properties of a simple two-dimensional model whose ground 
state is a d-wave superconductor. Using a self-consistent, conserving formulation, we show that 
pairing correlations above Tc lead to the appearance of a highly anisotropic pseudogap in the elec- 
tronic spectral function and the destruction of the Fermi surface. We discuss the similarities and 
differences between our results and ARPES experiments on underdoped cuprates. 

PACS numbers: 74.25.-q, 74.20.-z, 74.20.Mn, 74.25.Dw 



The deviations from Fermi liquid theory (FLT) in the 
normal state of high Tc superconductors are now well 
established. It has recently become clear that the under- 
doped cuprates exhibit even more remarkable deviations 
from FLT than the optimally doped materials: not only 
are the quasiparticles not defined, but the Fermi surface 
also becomes ill-defined due to the opening of a pseudo- 
gap. Early evidence for the suppression of low frequency 
spectral weight above Tc came from a variety of probes, 
including NMR , specific heat [|| and optics |^ . Recent 
angle-resolved photoemission (ARPES) studies [^Sj on 
Bi-2212 compounds have considerably clarified the situa- 
tion by providing direct evidence for a highly anisotropic 
pseudogap, for Tc < T < T* , which is similar in its 
magnitude and in its angular dependence to the d-wave 
superconducting gap below Tc. 

These observations of a normal state pseudogap find 
a simple explanation in a theory in which the pair- 
ing amplitude develops at a crossover scale T* higher 
than the Tc ^ ns/m* at which phase coherence sets 
in. The separation between Tc and T* naturally occurs 
in low density, short coherence length superconductors, 
and leads to striking deviations from FLT in degenerate 
Fermi systems in 2D Q]. However, rather little is known 
theoretically about the pseudogap state above Tc in d- 
wave superconductors (since, for technical reasons, the 
quantum Monte Carlo results are restricted to s-wave 
pairing). Such a study is clearly important, not only 
because the experiments show d-wave pairing, but also 
to compare the predictions of these theories |^,0J|] with 
those of RVB-based theories [|lO| in which spinous pair at 
T* and holons condense at Tc- An improved understand- 
ing of the origin of this pseudogap effect may also be 
an important clue towards a theory of high temperature 
superconductivity. 

As a first step in this direction, we study a phenomeno- 
logical model of a 2D d-wave superconductor, in a param- 
eter range where the normal state is dominated by d-wave 
pairing correlations. Using a self-consistent, conserving 
aproximation, we show that there is a crossover tempera- 



ture scale T* below which normal state spectral functions 
exhibit anomalous dispersion with a highly anisotropic 
suppression of spectral weight - pseudogap - near the 
chemical potential. This also leads to the partial destruc- 
tion of the Fermi surface along certain directions in the 
Brillouin zone. Both these effects are very similar to the 
ARPES experiments. We also find that the normal state 
spin susceptibility acquires a spin-gap like T-dependence. 

Let us consider a simple two-dimensional model which 
has a superconducting ground state with d^^-y^ symme- 
try, defined by 
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where the d-wave separable potential is 14:, k' — 
Udfi'k)f{k') with Ud < and /(k) = {coskx - cos fey). 
This potential is a piece of the nearest-neighbour interac- 
tion used in Ref. [ pd] | to model superconductivity in the 
cuprates. For simplicity, we study the nearest-neighbor 
dispersion Ck = —2t(cosk,j: + cosfcj,), with i = 1, on a 
square lattice with periodic boundary conditions. The 
chemical potential /i is adjusted to obtain the required 
density n. 

To investigate the finite temperature properties of this 
model in the intermediate coupling regime {\Ud\ of order 
bandwidth), we use the "fluctuation exchange" (FLEX) 
approximation [ p2[ . Since the important correlations in 
this problem are in the particle-particle (p-p) channel, 
we dress the propagator with these and solve the prob- 
lem self-consistently. Specifically, the vertex (with two 
incoming and two outgoing legs) is defined by the stan- 
dard integral equation written symbolically as fl^ : 
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where j., is the p-p irreducible vertex and / is the 
same quantity with twisted outgoing legs. Here and be- 
low all the quantities are matrices in spin space, the 
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four-momentum k = (k, ifc„) with ifc„ a fermion Mat- 
subara frequency {p and k' have the same nature) and 
the four-momentum q = (q, with ivrn a boson Mat- 
subara frequency, and the symbohc summation means 
integrating out (with proper factors) intermediate mo- 
menta and frequencies and matrix multiphcation with 
respect to spin indices. The Green's function G{k) = 
[ikn — Ck — M — 5](fc)]^^ is defined in terms of the self- 
energy, which satisfies the relation [0| 



I](fc) = J2 Vk,pG{q^k)G{p)G{q-p)Tl 
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where Vk^p is the bare two-body potential. 

We now approximate the p-p irreducible vertex by the 
bare potential from (]l|), in which case the contribution 
from / to S vanishes. For the self-energy, this is com- 
pletely equivalent to the well-known self-consistent T- 
matrix approximation [l4| . The advantage of the present 
approach is that in addition to the one-particle Green's 
function F also defines the two-particle Green's functions, 
and diagrams with "twisted" legs are important for eval- 
uating response functions such as the spin susceptibility 
(as we will discuss in detail elsewhere |15 ). 

The coupled integral equations for the self-energy and 
the vertex part are numerically solved using fast Fourier 
transforms [isl and the analytic continuation from Mat- 
subara to real frequencies is performed using Fade ap- 
proximants Technical details, as well as checks on 

the method and the numerics for the case of the attrac- 
tive Hubbard model, where we could compare our results 
against quantum Monte Carlo [Q, will be described else- 
where [ |l5| . 

In this paper we focus primarily on the spectral func- 
tion A(k,uj) = — (l/7r)ImG'(k, w -I- «0+), which is closely 
related to the ARFES intensity |18|. We show repre- 
sentative results for Ud = —8 and n = 0.5 (quarter- 
filling. Ha = 0.25); qualitatively similar results were 
found for several other parameter sets. All the results 
are in the non-superconducting state above Tc (see be- 
low) . In Fig. |l| we plot the spectral functions at a high 
temperature T = 2.0 for k varying from (0,0) to (tt, 0), 
and see the peaks of A(k, oj) disperse through the chem- 
ical potential {ui = 0). We may identify points k* in the 
Brillouin zone such that A(k*,w) has a dominant peak 
at w = 0. The "locus of gapless excitations" {k*} then 
generalizes the notion of a "Fermi surface" (FS) to 
finite temperatures without any assumptions about well- 
defined quasiparticles, and, quite generally, it is a closed 
contour in the repeated-zone scheme. 

It is worth commenting on the lineshapes in Fig. |l|: at 
(0, 0) the peak is actually infinitely sharp for our model, 
and the small width is put in by hand. Our choice of Vk_k' 
implies that states along the diagonal (0, 0) to (tt, tt) are 
totally unaffected by interactions. On the other hand, 
this potential gives rise to very strong effects near (tt, 0) 
where the spectral functions acquire very large widths 
as seen from Fig. |^, so that the quasi-particle nature 
is completely destroyed. Note that a log scale is used 
so that this spreading of spectral weight over an enor- 
mous frequency range can be easily seen. The sum rule 
J du!A{\i., cij) = 1 is satisfied (to very high precision) for 
each k. 
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FIG. 1. A(]i.,uj) for a sequence of momenta k = {xn/32, 0); 
X = {0, 6, 9, 12, 15, 18, 21, 24, 27, 32} for Ua = -8 and n = 0.5 
at a high temperature T = 2.0. In the inset the points in- 
dicate the momenta k and the solid curve shows the T = 0, 
non-interacting FS. 

As the temperature is lowered below a scale T* ~ 1 
(for the same choice of parameters as in Fig. |^), remark- 
able changes takes place in the spectra, as seen from 
Fig. H, which shows results at T = 0.2 (but still above 
Tc). First, we find that spectra which showed one broad 
feature at high T now show a multiple peak structure (see 
further below). Second, the dominant peaks of ^(k, w) 
exhibit very anomalous dispersion: as k varies from (0, 0) 
to (tt, 0), the peak approaches a; = but never crosses it, 
either "bouncing back" towards negative w or "stalling" 
(depending on parameter values), in complete contrast 
to the high T results described above. This is exactly 
like the pseudogap behavior seen in ARPES experiments 
on underdoped cuprates in the temperature regime 

Tr<T< T*. 
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FIG. 2. j4(k,aj) for a sequence of momenta k = {xn/32, 0); 
X = {0, 6, 9, 12, 15, 18, 21, 24, 27, 32} for Ua = -8 and n = 0.5 
at T = 0.2. 

To understand the multiple peak structure of the spec- 
tral functions it is useful to look at plots of the real and 
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imaginary parts of the (retarded) self-energy as functions 
of Lo. We find that the dominant peak of A{k,uj), which 
is at w < for each k in Fig. ||, is associated with a solu- 
tion of cj—ek—/i—ReE(k, w) = for which dReE/dio < 
and ImS is small. The very broad features at lu > come 
from solutions for which ReS has a positive slope and 
IniE is large. It is tempting to describe the "bounce" and 
multiple-peak structure as a precursor of the Bogoliubov- 
like dispersion of excitations in the SC state. In the pres- 
ence of strong self-energy effects, such a simple picture 
may need to be generalized significantly. 

Let us now ask how the pseudogap affects the "Fermi 
surface" (FS) by studying the "locus of gapless excita- 
tions" {k*} (defined, as before, by the condition that the 
spectral function at that k has a dominant peak cen- 
tered at a; = 0). Along (or near) the zone diagonal, 
(0,0) to (tt, tt), interaction effects are absent (or weak) 
and there is a conventional FS crossing with a well de- 
fined k*. However, along the (0,0) to (7r,0) there is no 
FS crossing. By studying the spectral function in the 
entire zone we find that the anomalous dispersion and 
large line-widths destroy the notion of a FS as a closed 
contour of gapless excitations (even in the repeated-zone 
scheme) . The resulting picture emerging from our calcu- 
lations, and consistent with ARPES experiments P,|5|,p0t , 
is that the Fermi Surface is destroyed in patches in the 
Brillouin zone, as schematically depicted in Fig. |^. 




FIG. 3. Fermi Surface destroyed in patches; solid lines rep- 
resent gapless excitations and the shaded patches indicate mo- 
menta where there is strong scattering and the FS is washed 
out. 

A more quantitative understanding of the destruction 
of the FS can be obtained by studying the angular depen- 
dence of the pseudogap and its variation with tempera- 
ture. In the spirit of the ARPES experiments we esti- 
mate the pseudogap by making scans through k-space 
and noting the position of the spectral function peak 
which is farthest to the right, i.e. at the largest frequency 
below zero. We then plot the spectral function pseudo- 
gap as a function of = arctan(fcj,//ea;). In Fig. W 
we plot the angle-dependence of Aps{9) at two tempera- 
tures T = 0.2t and T = 0.75t. The first point to notice 
is the strong anisotropy of the gap, which is always sup- 
pressed to zero in an arc about the diagonal. The extent 
of the zero gap (nodal) region , and the magnitude of 
the maximum gap at low T are both sensitive functions 
of the choice of parameters (see below); we find that the 
larger the maximum gap, the smaller the nodal region. 
The second important point to note is the T-dependence 
of the pseudogap, which suggests a gap collapse due to 
quasiparticles excited around the nodal regions in a d- 
wave SC gap, similar to a recent suggestion of Lee and 



Wen At high temperatures, the pseudogap gets sup- 
pressed and eventually disappears upon further heating; 
above a crossover scale T* 2± 1 (for the parameters dis- 
cussed here) the spectral function peaks disperse through 
uj = for all fixed-fcj^ scans. Thus we find that we re- 
cover a closed contour of gapless excitations above T*. 
However, in our model we find that this "Fermi surface" 
has T-dependence It is quite remarkable that in 

the ARPES experiments, there are indications ||^ of an 
underlying Luttinger FS which is T-independent within 
error bars. 
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FIG. 4. The angle dependence of the Aps{6) for tempera- 
tures T = 0.2 and T = 0.75 and Ud = -8, n = 0.5. 

It is instructive to consider the doping dependence of 
the effect. In Fig |^ we show Aps{9) at the same tem- 
perature and coupling, for three values of the density. 
The pseudogap clearly increases rapidly with increasing 
density. 
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FIG. 5. The angle dependence of Aps{0) for densities 
n = 0.5, n = 0.55 and n = 0.60 and Ud = -8, T = 0.2. 



Up to now, we have discussed the pseudogap behaviour 
in our model in the single-particle spectral weight. In 
Fig. ^ we show that the temperature dependence of the 
static, uniform, spin susceptibility also has a spin gap. 
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similar to what is observed in NMR and has been cal- 
culated [0 for the case of the attractive Hubbard model 
which has s-wave pairing. More details of this calculation 
will be presented elsewhere . 
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FIG. 6. The temperature dependence of the spin suscepti- 
bility for Ud = —8 and Ud ~ —4 at density n — 0.5. 

To ensure that we are above the superconducting Tc 
we calculate the parameter A \Ud\XppiQ = 0) where 
X^piq) = {T/N)Y.kP{k)G{k)G{q -k). The gap equa- 
tion at T — Tc reduces to the condition A = 1, which 
in our approximation coincides with the appearance of 
a pole in the full vertex. Thus, one should keep A < 1 
to be in the normal state. In 2D, there is of course no 
LRO at non-zero temperatures and one expects A < 1 for 
all T > 0. Even though our formalism does not include 
topological excitations, it is reasonable to expect Awl 
near the onset of algebraic order. In practice we consider 
only temperatures for which A < 0.8. 

We conclude this paper by showing that the pseudo- 
gap behavior described above is occuring in a degenerate 
Fermi system above Tc- This is important to establish 
since gap-like features can be trivially obtained either 
in a system below Tc, or in a strongly coupled regime 
where all the electrons are tightly bound up into bosonic 
pairs. The specific heat |^ and the ARPES |^,|| exper- 
iments clearly indicate that in the underdoped systems 
one is still dealing with a degenerate Fermi system and 
not bosons. By looking at the momentum distribution 
n(k) jl^], and by ascertaining that the chemical poten- 
tial fj, ^ T, we know that we are in a degenerate Fermi 
regime. 

In this paper we have studied the normal state of a 
simple model in which d-wave pairing correlations above 
Tc lead to the appearance of a highly anisotropic pseu- 
dogap and the destruction of the Fermi surface, which 
are remarkably similar to the results of ARPES exper- 
iments. Quantitative comparison with the experiments 
must, however, await a controlled calculation based on a 
microscopic model which describes how a Mott insulator 
upon doping goes into a short coherence d-wave super- 
conductor whose normal state is dominated by pairing 
correlations. 

A. N. and J. E. are supported by Boston College. Addi- 



tional support by the National High Magnetic Field Lab 
and Los-Alamos National Lab is acknowledged. J. E. 
and M. R. would like to thank Professor Yu Lu and the 
ICTP, Trieste, for hospitality during a visit which helped 
to initiate this research. 



[1] W. W. Warren et al, Phys. Rev. Lett. 62, 1193 (1989); 

M. Takigawa et al, Phys. Rev. B 43, 247 (1991). 
[2] J. Loram et al, Physica C 235-240, 134 (1994). 
[3] L. D. Rotter et al, Phys. Rev. Lett. 67, 2741 (1991); C. 

C. Homes et al, Phys. Rev. Lett. 71, 1645 (1993). 
[4] A. Loesser et al, Science, bf 273, 325 (1996); J. Harris et 

al, Phys. Rev. B 54, R15665 (1996). 
[5] H. Ding, et al. Nature 382, 51 (1996) and Phys. Rev. 

Lett. 78, 2628 (1997). 
[6] C. Sa de Melo, M. Randeria, and J. R. Engelbrecht, Phys. 

Rev. Lett. 71, 3202 (1993). 
[7] M. Randeria et al, Phys. Rev. Lett. 69, 2001 (1992); 

N. Trivedi and M. Randeria, Phys. Rev. Lett. 75, 312 

(1995). 

[8] Y. J. Uemura et al, Phys. Rev. Lett. 62, 2317 (1989). 
[9] V. Emery and S. Kivelson, Nature 374, 434 (1995). 

[10] G. Kotliar and J. Liu, Phys. Rev. B 38, 5142 (1988); H. 
Fukuyama, Prog. Theor. Phys. Suppl. 108, 287 (1992); P. 
A. Lee and X. G. Wen, Phys. Rev. Lett. 76, 503 (1996). 

[11] A. Nazarenko, A. Moreo, E. Dagotto and J. Riera, Phys. 
Rev. B 54, R768 (1996). 

[12] N. E. Bickers and D. J. Scalapino, Ann. Phys. (N.Y.) 
193, 206 (1989); N. F. Berk and J. R. Schrieffer, Phys. 
Rev. Lett. 17, 433 (1966). For applications to the repul- 
sive Hubbard model, see, e.g., J. W. Serene and D. W. 
Hess, Phys. Rev. B44, 3391 (1991) and C-H. Pao and N. 
E. Bickers, Phys. Rev. B51, 16310 (1995). 

[13] See, e.g., "Theory of Finite Fermi Systems and Applica- 
tions to Atomic Nuclei" by A. B. Migdal (Wiley Inter- 
science, NY, 1967). 

[14] R. Haussmass, Phys. Rev. B 49, 12975 (1994); R. Micnas 
et al, Phys. Rev. B 52, 16223 (1995). 

[15] A. Nazarenko, J. R. Engelbrecht, M. Randeria, and E. 
Dagotto, (in preparation). 

[16] The FFT works reliably with functions that decay as 
1/cj^ (or faster) at large ui because of the necessity of 
a cutoff. For the present model a cutoff 10 to 20 times 
the bandwidth was necessary. Our runs were typically for 
64x64 lattices with a uniform uu mesh of 512 points. 

[17] H. Vidberg and J. Serene, J. Low Temp. Phys. 29, 179 
(1977). 

[18] See: M. Randeria et al, Phys. Rev. Lett. 74, 4951 (1995) 
and references therein. 

[19] The gapless locus for a strongly correlated system at fi- 
nite temperatures, need not posses all of the properties 
that we associate with the usual T = definition of a 
FS in a Fermi liquid. For instance, there is no guarantee 
that the volume enclosed by this locus is T-independent. 
Nevertheless, by focusing on gapless excitations, it seems 
to be a very natural generalization. 

[20] Z.-X. Shen and J. R. Schrieffer, Phys. Rev. Lett. 78, 1771 
(1997). 

[21] P. A. Lee and X. G. Wen, cond-mat preprint (Feb. 97). 



4 



